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We report on advances in the study of the cores of NGC 6302 and NGC 6537 using 
infrared grating and echelle spectroscopy. In NGC 6302, emission lines from species 
spanning a large range of ionization potential, and in particular [Siix] 3.934 /im, are 
■ interpreted using photoionization models (including CLOUDY), which allow us to 

| reestimate the central star's temperature to be about 250000 K. All of the detected 

£SJ , lines are consistent with this value, except for [Al v] and [Alvi]. Aluminium is found to 

be depleted to one hundredth of the solar abundance, which provides further evidence 
for some dust being mixed with the highly ionized gas (with photons harder than 
ON . 154 eV). A similar depletion pattern is observed in NGC 6537. Echelle spectroscopy 

of IR coronal ions in NGC 6302 reveals a stratified structure in ionization potential, 
^ (— | which confirms photoionization to be the dominant ionization mechanism. The lines 

O are narrow (< 22 kins" 1 FWHM), with no evidence of the broad wings found in optical 

' \ lines from species with similar ionization potentials, such as [Ne v] 3426 A. We note 

the absence of a hot bubble, or a wind blown bipolar cavity filled with a hot plasma, 
at least on 1" and 10 kms -1 scales. The systemic heliocentric velocities for NGC 6302 
and NGC 6537, measured from the echelle spectra of IR recombination lines, are found 
to be -34.8±1 kms -1 and -17.8±3 kms -1 . We also provide accurate new wavelengths 
for several of the infrared coronal lines observed with the echelle. 

' Key words: planetary nebulae: NGC 6302: NGC 6537 - infrared: ISM: lines and 

^ , bands - ISM: abundances - ISM: kinematics and dynamics. 



1 INTRODUCTION 

The extreme high excitation PNe NGC 6302 and 
NGC 6537, which are both bipolar and of Peimbert's 
Type I, show lines from the highest ionization species 
found in PNe (e.g. [Sivi] 1.96 jum, Ashley & Hyland 
1988, hereafter AH88) as well as fast flows of order 
1000 kms -1 (Meaburn & Walsh 1980). Proposed ioniza- 
tion mechanisms include photoionization by a very hot 
central star, or shock excitation (e.g. Rowlands et al. 
1994, Lame & Ferland 1991, AH88, Barral et al. 1982). 
Pottasch et al. (1996) published a portion of the ISO 
[Infrared Space Observatory) spectrum of NGC 6302, 
and found that the properties of the spectrum can be 
well accounted for by photoionization by a 380000 K 
greybody central star, using the photoionization code 
CLOUDY (Ferland 1996). 

We report here on a study of the cores of NGC 6302 
and NGC 6537 using near- and mid-infrared spec- 



troscopy. Low and medium resolution CGS3 and CGS4 
infrared spectra taken at UKIRT have higher sensitivity 
in the atmospheric windows than the SWS on ISO, and 
our spectrum of NGC 6302 shows the presence of [Si ix] 
3.934/im, which had only been tentatively detected in 
the ISO spectrum. This detection has allowed us to 
reestimate the central star's temperature to be about 
250 000 K, the upper limit inferred by Marconi et al. 
(1996) from the non-detection of [Siix] by Oliva et al. 
(1996). 

We also present unprecedently high dispersion in- 
frared echelle spectra of the cores of these nebulae. 
Knowledge of the physical conditions and velocity fields 
in the obscured cores of bipolar PNe is the key to con- 
structing models, but these are currently inferred only 
through extrapolation from optical studies. We show 
here that tracing the velocity field with coronal ions 
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leads to quite a different picture of the structure of these 
bipolar planetary nebulae. 

The infrared coronal lines are in fact very nar- 
row, with no evidence of the fast flows detected in 
[Nev] 3426 A towards the core of NGC 6302 by 
Meaburn & Walsh (1980). We confirm the inference by 
AH88 that the full width at zero intensity of [Sivi] is 
less than ^30 km s -1 . The absence of broad wings in the 
infrared lines constitutes another piece of information 
with which to understand the physical conditions in the 
regions where the [Nev] 3426 A broad wings originate. 

The CGS4 spectra allow the determination of sev- 
eral coronal line wavelengths, which we measure to be 
significantly different from the ISO results (Feuchtgru- 
ber et al. 1997). We present a thorough account of our 
wavelength calibration procedures, and a discussion of 
their accuracy. 

We describe our observations in Section |[ the 
analysis of the spectra and the determination of ac- 
curate wavelengths for several coronal ions are de- 
ferred until the Appendix. Section ^ consists of the 
model-dependent interpretation of our observations. 
We present an abundance analysis for NGC 6302 and 
NGC 6537 and estimate the central stars' temperatures. 
In that Section we also use narrow band images of 
NGC 6302 to infer properties of the progenitor. Section 
^ contains an interpretation of the kinematical informa- 
tion from the echelle results, while Section || summarises 
our conclusions. 



2 OBSERVATIONS 

2.1 Low resolution spectroscopy 

The 8-13 fjxa spectra (atmospheric N-band window) 
of NGC 6302, PNG349. 5+01.0, and NGC 6537, 
PNG010.1+00.7 were obtained with CGS3 at UKIRT, 
on May 15 1994. The observations were chopped for 
background cancellation. The 32 Si: As detectors and the 
high resolution grating provide a resolution of 0.05/im 
(R=200 at 10/xm), sampled 3 times per resolution ele- 
ment. The reduced spectra, obtained through a 4" cir- 
cular aperture, are shown in Figure |l|. 

NGC 6302 and NGC 6537 were also observed at 
UKIRT through the atmospheric windows centred on 
the L and M filter band-passes, on April 24, 25 and 
26 1994, using CGS4 equipped with the short camera. 
The observations were carried out in chopping mode to 
allow background cancellation, at a resolving power of 
R=1500, a spatial resolution of 3" along the slit and 
a slit width of 3", with the spectra oversampled three 
times by shifting the detector's position 6 times over two 
pixels. The CGS4 detector at the time was a 62x58 InSb 
array. Due to problems with the array electronics pro- 
ducing random bad pixels in the first 20 columns of the 
array, the spectra were extracted using the minimum 
number of rows necessary to cover most of the emission 
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Figure 1: CGS3 spectra of NGC 6302 and NGC 6537. 
The spectrum with annotations is that of NGC 6302. 
Wavelengths are in microns, and flux densities are in 
10~ 12 W mr 2 plvo.- 1 . 

from the nebula. The reduced spectra are shown in Fig- 
ure |^. The dominant features are fine structure emis- 
sion lines from high stages of ionization, which probe 
the cores of the nebulae. 



2.2 Echelle spectroscopy. 

The nebulae were observed with the CGS4 short camera 
and echelle at UKIRT on May 11 and 12 1997. Two out 
of three nights were clear, so the observations are not as 
complete as hoped. In the case of NGC 6302, an infrared 
disk lies at the waist of the bipolar structure, where a 
dark lane is seen in optical images (Lester & Dinerstein 
1984). We oriented the CGS4 slit along two perpendicu- 
lar position angles, the || position lies at —20°, along the 
dark lane at the waist of the nebula, and the _L position 
lies at -110°, along the bipolar axis. We only obtained 
echelle spectra with the || position angle for NGC 6537, 
at -45°. The list of observed lines and their orientation 
in NGC 6302 is: Br 7 , _L & ||; [Mgvm], _L & ||; [Sivi], ||; 
[Arvi], _L. In NGC 6537, we observed Br7, [Sivi], and 
[Alv], all at the || slit orientation. 

At wavelengths near 4/j,m, the order sorting CVF 
filter produced interference fringes, hampering the data 
acquisition. In the case of [Arvi], the line was close to 
a fringe maximum, and the spatial distribution could 
still be determined. This ionic line was one of the best 
candidates to trace the velocity field since it is one of 
the very brightest, and also because its ionization po- 
tential is only a little less than that of Ne v, in which 
Meaburn & Walsh (1980) detected 1000 kms -1 wide 
wings. We were not able to obtain echelle measurements 
of [Siix] 3.934/mi, the highest excitation line so far de- 
tected from NGC 6302. Transmission modulation by the 
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Figure 2: CGS4 grating spectra of NGC 6302 (upper) and NGC 6537 (lower). Wavelengths are in microns, and 
flux densities in 10 -12 W m -2 /im" 1 . Superposed on NGC 6537's spectrum, in light grey, is a low resolution model 
atmospheric transmission function appropriate for the summit of Mauna Kea. 



fringes completely dominates the emission from [Si ix] , 
which is more than a thousand times fainter than the 
[Arvi] line. 

A resolving power of ^20 000 can be achieved with 
the echelle and UKIRT's CGS4 spectrometer, corre- 
sponding to a velocity resolution of 15 km s _1 . With 
the short camera and a one pixel wide slit, the pixel 
size is about 1.5" spatially and 1" spectrally. We over- 
sampled the spectra by shifting the array 3 times over 
2 pixels. 

Figure |^ and |^ show the spectra obtained by col- 
lapsing all the emission along the slit. In the cases where 
it was relevant (i.e. when there was significant struc- 
ture), the spectrum of a standard star is superposed in 
dotted lines on the nebular ion's spectrum. The stan- 
dard star spectra are shown in arbitrary units, but zero 
counts is represented by the horizontal dotted lines. Fig- 
ures ^| and [l^ contain the reduced frames for the lines 
observed with CGS4+echelle. 

The secondary line at 3.0271/im apparent in the 
[Mgvm] spectra for both slit positions could not be 
identified. It is not believed to be a [Mgvm] 3.028/im 
velocity component at —64 kins" 1 because it is not seen 
in the spectra of the other ions. Neither is it an atmo- 
spheric transmission feature, which can be concluded 
by inspecting the frames presented in Figure ^|, nor 
is it a leak from neighbouring orders. It is nonethe- 
less probable that the secondary line originates from 
a coronal ion formed with an ionization potential be- 
tween that required to produce Si VI and Siix (i.e. from 
166.77 eV to 303.16 eV), since the low resolution spec- 
trum of NGC 6537 shows no sign of any line at all 



near 3.028^m. A Doppler shift of the secondary line 
deeper into the atmospheric transmission trough seen 
in Figure || is discarded on the basis that at the time of 
the echelle observations, the velocity of NGC 6302 was 
— 49±1 kms" 1 , in the Earth's rest frame, whereas at 
the time of the low resolution observations, NGC 6537 
was at — 42±3 kms -1 (see Appendix |X|). 

In the Br7 frame (Figure |3|) , the fluxes and wave- 
lengths of the two secondary lines correspond to He I 4- 
7 at 2.16415^m (the singlet 1 F°- 1 G and triplet 3 F°- 3 G 
transitions are blended) and Hell 8-14 at 2.16464/^m, 
using the emissivity Tables from Smits et al. (1991) and 
Hummer & Storey (1987). 

Only the raw observation is shown for [Ar vi] . This 
frame was affected by the CVF fringes, which can be 
seen in the standard star's spectrum shown by the dot- 
ted line. However a flux calibration is possible and the 
total [Arvi] flux through the one-pixel 1" wide slit 
is about 10 times that for [Mgvm]. A higher [Arvi] 
flux relative to [Mgvm] is found in the low resolution 
spectrum, with a 3" slit, which is an indication that 
the [Arvi] emission has a larger spatial extent than 
[Mgvm]. 

3 ABUNDANCE STUDY AND 
PHOTOIONIZATION MODELS 

3.1 Extinction estimates and emission 
measures 

We shall adopt for NGC 6302 the density and temper- 
ature estimated by Oliva et al. (1996) from a variety of 




Figure 3: Lines observed from NGC 6302 with the CGS4 echelle. Flux densities are given in 10 -12 Wm" 2 /an, 
and the calibrated wavelengths for the nebula's rest frame are in microns. The dotted lines show the standard star 
spectra used to flux calibrate the observations, doppler shifted to the rest frame of NGC 6302. The feature in the 
short wavelength wing of the [Mg vni] line could not be identified. The dotted curve over the [Ar vi] spectrum 
corresponds to the CVF fringe pattern. The rms noise level is 5 x 10 -15 Wm -2 fim in the [Mg vin] spectra, and 19 
counts in the [Arvi] spectrum. 
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Figure 4: Lines observed from NGC 6537. Flux densities 
are given in 10 -12 W m~ 2 /im, and wavelengths are in 
microns and in the nebula's rest frame. The dotted lines 
are the standard star spectra used to flux calibrate the 
observations, doppler shifted to NGC 6537's rest frame. 



transitions, T e =19 000K (e.g. from [Arv] A4626,A7006 
and [Neiv] A2422,A4725), and AT c =18000cm~ 3 . For 
NGC 6537 we used the [O in] electron temperature 
T e =16000K and A o =17000cm~ 3 from Feibelman et 
al. (1985). In the case of a photoionization dominated 
structure the temperature profile across a uniform neb- 
ula is expected to be fairly flat around T e =20000 K 
(CLOUDY models at constant proton density, see be- 
low) . We will also argue in Section Q in a comparison 
between the spatio-kinematical data and the structure 
of model nebulae, that the physical conditions must be 
quite uniform across the cores of the nebulae, at least 
when averaged over 1" scales. 

We used the relative intensities between H I recom- 
bination lines to evaluate errors in flux calibration. The 
errors stem from variations in slit positions for differ- 
ent frames, approximate estimations of standard star 
fluxes, and unresolved atmospheric transmission fea- 
tures. A restricted range in the J window including Pa/3 
at 1.282/im and the 3 3 S-4 3 P Hei transition at 1.2792 
was also observed in both nebulae with CGS4 in grat- 
ing mode. Atmospheric transmission is quite good for 
Pa/3, and it is a bright line observed at very good signal 
to noise. Pa/3 was therefore used as the flux reference 
H I line. On the other hand, Pa/3 is appreciably more af- 
fected by extinction than the longer wavelength IR lines, 
but its dereddened flux can be estimated from free-free 
radio continuum maps. We used the relation given by 
Milne & Aller (1975) to estimate the total H/3 flux from 
the 5 GHz free- free continuum measurement. We follow 
similar assumptions as in AH88, that the He abundance 
relative to Hi is 18 per cent, and the fraction of the He 
atoms that are doubly ionized is 75 per cent. The total 
Pa/3 flux can then be estimated from the emissivity ta- 
bles given by Hummer & Storey (1987), assuming Case 
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B. For the hydrogen lines with Hei blending we used 
the He I emissivities of Smits (1991), normalised to Hei 
3 3 S-4 3 P. 

Gomez et al. (1989) present a VLA map of 
NGC 6302 at 5 GHz, which allows extinction estimates 
for a 3x3 arcsec 2 beam. Our measured flux for Pa/3 of 
3.9xl0 -15 Wm -2 gives Ay = 3.7 in the central region, 
with S(5 GHz)=375mJy (over the same region we have 
1.6xlCr 16 Wm -2 for the Hei 3 3 S-4 3 P observed flux). 
Although this value is not in agreement with the mean 
value for the whole nebula of Ay = 2.65 (Pottasch et 
al. 1996), it agrees with the conclusion by AH88 that 
extinction in the central regions of NGC 6302 is higher 
than the mean by 1.3 magnitudes. Unfortunately no 
VLA maps are available for NGC 6537, but the esti- 
mations of AH88 give Ay = 4, very close to the value 
for NGC 6302 (for which they estimate Ay = 4.3). We 
fixed the extinction in NGC 6537 to Ay = 3.5 magni- 
tudes by requiring the unblended Bra to Pa/3 flux ratio 
to match the ratio of their emissivities at 20 000 K (from 
Hummer & Storey 1987). 

Unresolved atmospheric transmission (AT) troughs 
can be recognised from an atmospheric transmission 
model for Mauna Kea (based on the HITRAN database, 
Rothman et al. 1992). We took Pf5-10 as typical, since 
it lies in a trough, but its predicted flux was very close 
to the observed one. Unresolved atmospheric transmis- 
sion features are however a source of uncertainty in the 
low resolution flux measurements, as shown by AH88 in 
the case of [Si Vi] , which can be completely absorbed by 
atmospheric H2O. We used the extinction law given by 
Rieke & Lebofsky (1985), and the dereddened and de- 
blended H 1 fluxes are at most 20 per cent different from 
the values expected from the radio map in the case of 
NGC 6302. Observing from Mauna Kea, with a low col- 
umn of precipitable water vapour, minimizes the effect 
of AT troughs. 

The emission measure integrated along the line of 
sight, EM = f ds N p N c , can be estimated from the 
intensity of hydrogen emission lines (e.g. Osterbrock 
1989), 

I = J ds N p N e (hv a off )/(47r), (1) 

where a c ff is the effective recombination coefficient 
for Pa/3 (Hummer and Storey 1987). The Pa/3 
line flux from the central 3"x3" in NGC 6302 is 
F Pa/3 =3.9xlO~ 15 Wm~ 2 , giving EM = 5.6xl0 35 m -5 . 
For NGC 6537, F Pa/3 =3.8x 10~ 15 Wm" 2 , and EM = 
5xl0 35 m~ 5 . 

3.2 Abundance patterns and photoionization 
models 

Table |l| lists the deblended and dereddened surface 
brightnesses through the central CGS4 row (3" x 3" 
square aperture), together with the abundances calcu- 



lated using the atomic data referenced there. The [Si vi] 
and [Si vn] intensities were derived from UKIRT obser- 
vations on April 1986 with the CVF on UKT9, using a 
19.6 arcsec circular aperture. We also used some fluxes 
from Pottasch et al. (1996) and Rowlands et al. (1994), 
as indicated in Table 111 We estimated the average in- 
tensity in the central 3" x 3" by using the spatial in- 
formation from our grating and echelle CGS4 spectra 
(see the end of Appendix |A|) . The fluxes obtained with 
large beams were scaled by a factor 0.4 for NGC 6302 
and 0.58 for NGC 6537, to estimate the corresponding 
fluxes in the central 3" x 3" (except for [S ill] , for which 
we used a factor of 0.15 for NGC 6302, and 0.37 for 
NGC 6537). 

Although we list only ionic abundances, we can 
draw some conclusions about abundance patterns from 
those elements having several observed ionization stages 
and reliable line fluxes. Magnesium seems to be slightly 
depleted relative to solar, but not as strongly as alu- 
minium. Potassium shows a very small ionization frac- 
tion in K vii, when compared to Ne v at a similar ioniza- 
tion potential, or to Navi and NaVll. Although a some- 
what higher abundance in neon is expected in Type I 
PNe (Corradi & Schwarz 1993), it could be that potas- 
sium is slightly depleted, since no special overabundance 
is expected in sodium. 

We also used CLOUDY 90.03 (Ferland 1996) to con- 
struct a photoionization model of NGC 6302. The neb- 
ula is simulated as a spherically symmetric shell, al- 
though the fluxes of emission lines are calculated by 
truncating the shell in order to simulate the long slit ob- 
servations. The central star's radiation field is assumed 
to be that of a blackbody. No systematic depletion of 
metals on grains was assumed, principally because lit- 
tle is known about dust survival in PNe. Nonetheless, 
the abundances of elements likely to be depleted were 
allowed to vary, in particular magnesium shows an im- 
portant line that should be reproduced accurately. The 
abundances were otherwise taken as solar. We also var- 
ied the filling factor, the blackbody's temperature and 
the inner radius of the nebula. CLOUDY's optimisation 
driver was constrained by the fluxes of IR coronal ions 
and the H11 column density (estimated from the ob- 
served emission measure, with N e = 1.8x10 cm -3 ). 

The main problem in using CLOUDY for this ap- 
plication is that we are probing the central region of 
NGC 6302, whereas CLOUDY calculates the fluxes from 
the whole nebula. The cumulative fluxes versus radius 
could be matched with our spectra, leaving the cut- 
off radius as yet another free parameter. But one way 
around this problem is to reproduce the fluxes rela- 
tive to another ionic line. We took [Arvi] at 4.53 /j,m 
as the reference line, and the CLOUDY output can be 
seen in Table ||, as well as the best fit values for the 
free parameters. A total hydrogen volume density of 
log(./VH)=4.2 was given as input, together with a rela- 
tive helium abundance by number of 0.14 (using 0.18, as 
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Table 1: Ionic abundances by number for NGC 6302 and NGC 6537, relative to solar values. 



NGC 6302 NGC 6537 







A[/im] 


Idered 

Wm~ 2 sr 


— 1 


abundance d 


Idered 

Wm~ 2 sr 


— 1 


abundance d 


IP 

eV 


atomic 
data c 


Hi 


3-5 


1.282 


4.15x10" 






4 73 x 1 0" 


■5 




1 3 6 


I 


He i 


o3q_ 4 3p 


1 253 


5 89x10" 


-7 


47 


3 93x 1 0" 


■7 


9 68 


24.6 


2 


Hen 


6-7 


3 0Q1 7 


1.49x10" 


■5 


1 38 


1 62x10" 


•5 


1 60 


54.4 


I 


IV P TT 
IN C 11 


2p, ,„ 2p 
r 3/2 r l/2 




1 91 y1 0" 

_1 . Zi ± A 1U 


A 


n ^a 


A 94 v 1 n~ 


•5 


n 1 7 




q 



lMVi TTT 


3"p 3"p 


10.0 


7 \y 1 n- 

/./oXlU 


-4 a' 


1 

1 .z 








a 1 n 

41 .U 


4 


AT T T T 

No III 


3n 3y-> 


36.0 


4.27x 10 


5 a' 


0.9 








41.0 


4 


AT 

Ne v 


3n 3y-> 
Pi" "2 


14.32 


l. 42x10 


3 a 


0.36 


4.75x 10" 


4 fc 


0.13 


97.1 


5 


Ne v 


Per Pi 


24.31 


5.31x10" 


-4 a 


0.37 


O AA , , 1 A- 

2.60x 10 


4 fc 


0.19 


97.1 


5 


Nc vi 


2r> 2r> 
r l/2" r 3/2 


7.64 


1 a f\ « 1 1 n — 

l. 49x10 


-3 


0.18 


A OA, , 1 A- 

4.89x 10 


-4 


0.12 


126 


6 


Navi 


3r> 3n 
Pi- P2 


8.64 


O OA , , 1 A — 

2.20x 10 


-5 


0.33 


a r ^ , 1 a 

8.95x10 


6 


0.15 


138 


5 


Na VII 


2t> 2t> 
J^l/2" ^3/2 


4.690 


1 OA , , 1 A- 

1.89x10 


-5 


0.12 








172 


7 


TV T 

Mgv 


3 P2- 3 Pl 


5.608 


ry AA , , 1 A 

7.99x10 


-5 a 


0.11 








109 


4 


Mg vii 


3n 3y-> 
Pi" P2 


5.502 


6.15x10" 




A A , , 1 A — 2 

4.9x10 








187 


5 


Mg vin 


2t) 2n 

"l/2- "3/2 


3.028 


O rr ,, 1 a- 

8.55x10 


-6 


2.8x10 








225 


/? 




AI V 


2t3 2d 
r 3/2" *l/2 


2.883 


5.86X 10 


-8 


1 A w 1 A — 3 

1.0x10 








1 OA 

120 


3 


Al vi 


3 P,- 3 Pt 
a 2 ^ 1 


3.659 


6 66x10" 


■7 


1.9xl0~ 3 








154 


4 


Si vi 


2 P , 2 P , 
r 3/2" *l/2 


1.96 


1.11x10" 


4 


0.22 


5.88x10- 


■5 


0.13 


167 


3 


Si vn 


3 P 2 " 3 Pl 


2.48 


9.70x10- 


■5 


0.14 


1.59x10" 


-6 


2.6xl0~ 3 


205 


4 


Si ix 


3 Po- 3 Pi 


3.934 


1.42x10" 


-7 


3.0xl0~ 5 








303 


5 


Sm 


3 Pl" 3 P 2 2 


18.71 


4.17x10" 


■5 6 


0.11 


7.68x10- 


■5 b 


0.22 


23.3 


8 


Siv 


Pl/2" P3/2 


10.50 


2.08x10" 


4 


0.13 


2.91x10" 


-4 


0.38 


34.8 


9 


Arm 


3 P 2 - 3 Pl 


9.01 


7.66x10- 


■5 


0.38 


3.44x10- 


■5 


0.19 


27.6 


8 


Arv 


3 Pl- 3 P 2 


7.91 


4.29x10- 


■5 


8.4xl0~ 2 


2.91x10" 


■5 


6.2xl0~ 2 


59.8 


8 


Arv 


3 Po- 3 Pi 


13.11 


2.09x10- 


■5 


5.1xl0~ 2 


2.22x10- 


■5 


5.8xl0- 2 


59.8 


8 


Ar vi 


2 P , 2 P , 
r l/2" r 3/2 


4.528 


1.37x10" 


-4 


0.11 


3.04x10- 


-4 


0.32 


75.0 


10 


Km 


2 P 3 /2- 2 Pl/2 


4.616 


1.52x10" 


■6 


0.20 


5.09x10- 


-7 


7.3xl0- 2 


31.6 


11 


Kvn 


2 P , 2 P , 
r l/2" r 3/2 


3.189 


1.32x10" 


■6 


3.3xl0- 2 


9.83x10- 


-7 


2.7x10-2 


99.9 


10 



a ISO measurement, Pottasch et al. (1996); a ' Pottasch & Beintema (1997) 
b Rowlands et al. (1994). 

c The collision strengths were for the most part taken from the IRON project. 1, Hummer & Storey (1987); 2, Smits 
(1991); 3, Saraph & Tully 1994; 4, Butler & Zeipen (1994); 5, Lcnnon &Burke (1994); 6, Zhang (1994); 7, Ferland 
(1996, as used in CLOUDY); 8, Galavfs et al. (1995); 9, Johnson et al. (1986); 10, Saraph & Storey (1996); 11, Pelan 
& Berrington (1995) 

d Ionic abundance relative to the total solar abundance of each element, or (Nx/Np)/solar. Solar abundances from 
Grevesse and Anders (1989), Grevesse and Noels (1993). 
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Table 2: CLOUDY output from the best fit models for 
NGC 6302 and NGC 6537. Columns 2-5 list ionic line 
fluxes relative to [Arvi] 4.53/im = f.00. 



a) observed 



b) CLOUDY output 



filling factor 
inner radius 
blackbody T* 
luminosity 
Mg abundance 
Al abundance 



NGC 6302 
0.70 

6.2xl0 15 cm 
241 000 K 
8500 L Q 
f» 0.5 solar 
« 0.01 solar 



NGC 6537 
0.58 

3.8xl0 15 cm 
156 000 K 
8500 L 



ID X(fim) 


NGC 6302 
Model Obs. 




NGC 6537 
Model Obs. 




[Nev] 14.3 


5.7 


10 




3.6 


1.6 




[Nev] 24.2 


1.8 


3.9 




1.2 


0.86 




[Nevi] 7.64 


8.6 


11 




1.3 


1.6 




[Navn] 4.69 


7.6 10" 2 


0.14 










[Mgv] 5.60 


0.30 


0.60 










[Mgvn] 5.51 


0.31 


0.40 










[Mgvm] 3.03 


5.0 10" 2 


6.010" 


2 








[AlV] 2.9 


6.0 10" 4 


5.010" 


4 








[Alvi] 3.66 


4.8 10" 3 


5.010" 


3 








[Si vi] 1.96 


0.80 


0.80 




9.0 10" 2 


0.19 




[Sivn] 2.47 


0.32 


0.71 




2.5 10" 3 


5.210" 


3 


[Si ix] 3.93 


1.210" 3 


1.010" 


3 








[S iv] 10.5 


2.26 


1.52 




2.89 


0.96 




[Arv] 7.91 


9.3 10" 2 


0.30 




0.12 


9.610" 


2 


[Arv] 13.1 


7.6 10" 2 


0.15 




9.5 10" 2 


7.310" 


2 


[Ar Vi] 4.53 


1.00 


1.00 




1.00 


1.00 




[Navi] 8.64 


8.3 10" 2 


0.16 




1.7 10" 2 


2.910" 


2 


[K VH] 3.19 


5.2 10" 2 


1.010" 


2 


2.5 10" 2 


3.210" 


3 


N(Hn) a 


1.49 10 21 


1.5510 


n 


1.33 10 21 


1.4010 


21 



a Column density of H + in cm 2 integrated outwards from the 
central star. 



reported by Aller et al. 1981 and as derived in this work, 
does not change the results) . The total luminosity of the 
central star was fixed to be 8500 L©, simply because it 
is close to a standard PN value, although Pottasch et al. 
(f996) estimate L = 11000 L© for a distance of 1.6 kpc 
(Gomez et al. 1993). Using 11000 Lq leaves the model 
unchanged. The calculation stopped at an outer radius 
of 1.2xl0 17 cm because the package's lower limit elec- 
tron temperature of 4xl0 3 K was reached. 

The best fit black body temperature is about 
250000 K. This is to be compared to the very high 
values of 450000 K found by Lame & Ferland (1991), 
and 380 000 K by Pottasch et al. (1996). One important 
new datum used here is the [Si ix] flux, which is actu- 
ally quite low and would be significantly overestimated 
by models with T* higher than 250000 K. This result 
has already been mentioned by Marconi et al. (1996) in 
their study of the Seyfert 2 galaxy NGC 1068. Arguing 
against the production of Si IX by hot stars, they take as 
an example the nucleus of NGC 6302, and model the ab- 
sence of [Si ix] from their NGC 6302 spectra by putting 
an upper limit on of 250 000 K (using CLOUDY as 
well). The central star in NGC 6302 is thus probably 
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Figure 5: Ionization distribution for NGC 6302. 
The abscissa is the photon energy required to ion- 
ize the lower stage of ionization, and the ordinate is 
N(X +i+1 )/N(X +i ). The uncertainty in the observed ra- 
tios can be up to a factor of 2. 

not as hot as previously thought, T^=380 000 K overes- 
timates the [Si ix] flux by a factor 60 (with an optimised 
filling factor of 0.16, and an inner radius of 3.1 xlO 16 
cm). But could silicon be depleted? The CLOUDY model 
shows that 70% of all silicon nuclei are distributed in 
Si v, Si vi and Si VII. The match to the [Si vi] and [Si vn] 
lines should provide a good grip on the silicon abun- 
dance, and they are approximately reproduced using a 
solar silicon abundance. 

Together with the model results for NGC 6302, Ta- 
ble |^ shows the result of a CLOUDY model that re- 
produces the lines observed from NGC 6537, with a 
central star temperature of 156 000 K (using similar in- 
puts as for NGC 6302, with log(A H ) = 4.2; the calcula- 
tion stopped because the package's lower limit electron 
temperature was reached). [Alvi] was not detected to- 
wards NGC 6537, but its predicted flux assuming a so- 
lar Al abundance would be 5.7xl0 -2 times the [ArVl] 
flux (which is 6.4xl0" 14 Wm" 2 for the central 3" x 3" 
in NGC 6537). As the [Alvi] frame shows Hul7 close 
to the detection limit, with a flux 5.4xl0~ 4 times the 
[Arvi] flux, we can conclude than Al is depleted to be 
less than one hundredth solar. 

The abundance ratios of consecutive stages of ion- 
ization as a function of ionization potential are, in the 
case of photoionization, a probe of the high energy flux 
from the central star. Figure ||a shows the ionization dis- 
tribution calculated from the observed fluxes. We can 
use the CLOUDY model for NGC 6302 to construct a 
similar curve, shown in Figure |b. The curve is m re- 
markable agreement with the observations, except for 
the silicon ions. The Mg vin to Mg vn ratio also argues 
for a of about 250 000 K. Such good agreement is 
surprising considering all the approximations involved 
in the modelling. In a complete model of the nebula, 
a number of arbitrary assumptions about the recombi- 
nation data have to be made from the very first ioniza- 
tion stage. The abundance ratio of consecutive stages of 
ionization is inversely proportional to the total recom- 
bination coefficient, and little information is available 
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on the dielectronic part, especially at nebular temper- 
atures. But even if the atomic parameters were known 
with sufficient accuracy, as well as the temperature of 
the coronal line region, the central star's spectrum may 
differ significantly from a blackbody. 

With the cautions noted above (that a match be- 
tween model and observed fluxes is acceptable within 
a factor of <~2), the abundances computed by CLOUDY 
confirm the patterns from the direct estimate. For both 
NGC 6302 and NGC 6537, aluminium is depleted to one 
hundredth of the solar value, or even less. The magne- 
sium abundance is about half solar in NGC 6302. Neon 
seems to be overabundant in NGC 6302, but by less 
than a factor of two relative to the solar value, ft could 
be that potassium is depleted to one fifth of the so- 
lar abundance in NGC 6302 as well as in NGC 6537, 
although only one line is available. The uncertainties 
in the photoionization model may be reflected in the 
discrepant Si vn to Si vi abundance ratio (which is sen- 
sitive to extinction) and the model deficiency in [ArV] 
and in the sodium lines, but in all cases the discrepan- 
cies are less than a factor of 2-3. 

Allcr ct al. (1981) found that Ca was depleted to 
one tenth solar in NGC 6302, while Oliva et al. (1996) 
estimated Ca to be depleted by a factor ~20, and Fe 
by a factor of <~ 10 3 (see also Shields 1983 for deple- 
tion patterns in other PNe). It could be that titanium 
is also strongly depleted. According to the CLOUDY 
model, the flux of [Tivi] 1.7^m should be a little less 
than that of Br7, at about one fiftieth of [Arvi], but 
no spectrum in the H window has been published for 
NGC 6302. However, most of the titanium should be in 
Tivn or Ti VIII, and the ground state, magnetic dipole, 
[Tivn] 4144.2 A line is missing in the spectrum pub- 
lished by Allcr ct al. (1981). So it seems likely that Ti is 
strongly depleted. The strong depletion of Al, and pos- 
sibly of Ti, implies that some dust is well mixed with the 
gas phase, in the form of Ti and Al compounds, such as 
perovskite, CaTiOs, and corundum, AI2O3, which have 
very high sublimation temperatures. Gail and Sedlmayr 
(1998) argue that the first condensates in circumstellar 
outflows are Ti compounds, which could be nuclcation 
sites for the aluminium compounds and silicates. When 
the AGB ejecta are heated to nebular temperatures and 
exposed to hard UV radiation, grain ablation could ex- 
plain the observed depletions. Kingdon & Ferland (1997 
and references therein) suggest that strong ionic deple- 
tions are indicative of dust survival in the ionized phase 
of PNe. Examples arc Fc depletions inferred from Fe VII, 
which requires radiation harder than 100 eV. We note 
here that since 154 eV is required to ionize Alv, the 
Al depletion derived for NGC 6302 from Al VI provides 
evidence for dust survival in very hard radiation field 
environments. Pwa et al. (1984, 1986) also measured Al 
depletions of a factor ^100 for NGC 6543 and a factor 
-10 for BD+30°3639, from IUE observations of Aln 
and Aim absorption lines. 



It would seem that photoionization models can ac- 
count for the observed ionization distribution in broad 
terms. Zhekov & Pcrinotto (1996) attempted to model 
the IR coronal line emission from high excitation PNe 
for the case of a collisionally excited plasma. They 
used an interacting wind model, following Weaver et 
al. (1977), in which the slow AGB wind is accelerated 
into a shock by the fast white dwarf wind. The hot bub- 
ble grows as the shock expands and 'vaporizes' material 
from the cold shell. After working out the temperature 
structure of the bubble, they calculated its coronal line 
emission using the hot plasma emissivities. But their 
predicted fluxes arc in strong disagreement with those 
observed from NGC 6302, e.g. the predicted [Siix] and 
[Mgvm] fluxes should be of the same order according 
to their model. 

3.3 NGC 6302's progenitor 

In this Section we present arguments for NGC 6302 's 
progenitor having a lower ZAMS mass than previously 
thought. The total luminosity of the central star can 
be estimated from a sum of all flux measurements of 
NGC 6302 from the UV to the infrared (Pottasch et al. 
1996 and references therein). For a distance of 1.6 kpc 
(as determined by Gomez et al. 1993 using VLA ob- 
servations of the expansion rate), the central star's 
luminosity is about 1.1 xlO 4 L Q . We can now locate 
the star in the H-R diagram and compare with the 
theoretical tracks of Blocker (1995). A temperature of 
380000 K would locate the star on the track of a hy- 
drogen burning model with a 7M Q initial mass (and 
a core mass of Mh = 0.940 M Q ) 60 years or so af- 
ter leaving the AGB. But our estimate of 250 000 K 
would fall on a track intermediate between the Blocker 
(1995) models with (M ZA ms,M h ) = (4 M Q , 0.696 M Q ) 
and (5 M Q , 0.836 M ), or a 4-5 M Q initial mass, for ei- 
ther the hydrogen burning models or the born-again he- 
lium burners. It is difficult to distinguish between the 
two burning modes, as the time-scale is slower in the 
helium burners by a factor of 3 only. If the bulk of the 
nebular material was ejected at the time the star left the 
AGB, the age of about 2000 years necessary to account 
for NGC 6302's size of 10 17 cm at an expansion speed 
of 10 kms -1 may be better explained in terms of the 
born-again models of Blocker (1995). The PN progeni- 
tor could have had a ZAMS mass of ~4 M and have 
left the AGB on the verge of a thermal pulse, since the 
time-scale for a straight evolution from the AGB as a 
hydrogen burner is only 700 yrs (alternatively the star 
could have left the AGB just after a thermal pulse, and 
evolved on a similar time-scale as the born-again case, 
although Blocker does not provide tracks for this case). 

The OH maser observed in NGC 6302 by Payne et 
al. (1988) provides evidence for an oxygen rich environ- 
ment, as does the ISO detection of mid-IR crystalline 
silicate emission bands (Waters et al. 1996), and yet the 
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Figure 6: IRCAM image of the 3.3/im UIR-band emis- 
sion in NGC 6302. Pixel numbers are plotted along the 
axes. The pixel size was 0.6" and the width of the nar- 
row band filter was 0.3yLtm. The contour levels are .05, 
.1, .2, .3, .4, .5, .6, .7, .8, .9 of the peak flux, to allow 
comparison with the VLA 6 cm map from Gomez et al. 
(1989). North is down and East is to the right. 

nebula shows the family of 'unidentified infrared' (UIR) 
emission bands, attributed to carbon based molecules 
(Duley & Williams 1981, Leger & Puget 1984). The 
OH maser emission seems to come from the material 
surrounding the core of the nebula, since the cuts in ve- 
locity published by Payne et al. (1988) peak at different 
positions around the centre, over about 8". On the other 
hand the 3.3 fim UIR emission band is centrally peaked, 
as can be seen on our IRCAM frame shown in Figure 
|(| The extent of the UIR emission can be compared to 
the VLA map of Gomez et al. (1989). The UIR bands 
are typical of C-rich PNe (Cohen et al. 1986), and these 
observations support a scenario in which NGC 6302 's 
progenitor was an M giant that went through a carbon 
transition at the very end of its evolution, as is analysed 
in detail by Frost et al. (1998). However, maser amplifi- 
cation requires particular conditions, and therefore does 
not systematically trace the distribution of O-rich mate- 
rial, and the evidence for a stratification in C/O chem- 
ical balance is only marginal until high resolution OH 
observations are available and a model is constructed. 
Nonetheless NGC 6302 clearly includes both C- and O- 
rich environments. Values of C/H = 4.4xl0~ 4 , N/H = 
9.2xl0~ 4 and O/H = 5.0xl0~ 4 have been derived for 
NGC 6302 by M. J. Barlow, using the UV and optical 
relative line intensities measured by Barral et al. (1982) 
and Aller et al. (1981), respectively. Collisionally ex- 
cited lines were used to derive the abundances of the 



singly, doubly and triply ionized species, while the 6h- 
5g recombination lines in the 4600 A region were used to 
derive the abundances of the four-times ionized species, 
as in Clegg et al. (1987). 

Synthetic AGB evolution models (Renzini & Voli 
1981) are compilations of model results from full evolu- 
tionary calculations, and allow searches over a continu- 
ous range of progenitor masses. Using an adaptation of 
the synthetic AGB model by Groenewegen & deJong 
(1993), we can search for a progenitor mass that ac- 
counts for the properties of NGC 6302. In this model 
we do not constrain envelope burning to stop at the 
same time as dredge- up (e.g. Frost et al. 1998). Many 
uncertainties affect the treatment of hot bottom burn- 
ing (HBB), in this case one of them being the arbitrary 
onset of HBB for core masses greater than 0.8 Mq^| 
But our purpose here is to illustrate an alternative ar- 
gument to the post- AGB tracks for NGC 6302 having 
a relatively low progenitor mass. The above uncertain- 
ties in the treatment of HBB would decrease the lower 
mass limit of the range of initial masses that account 
for NGC 6302's composition. Figure ^ shows the time 
after the C transition as a function of initial mass for 
stars at the end of the AGB, as well as the averaged 
C/O and N/O ratios over the last 10000 years of mass 
loss. Towards 6M© the time after the C-transition is 
too long, the carbon rich nebular material having been 
ejected at most over 2000 years. Below 4.2 M Q HBB 
is not activated and the envelope composition is domi- 
nated by the dredge-up of carbon. Thus, in terms of the 
quenching of HBB at the end of the AGB evolution of 
the progenitor, the main sequence mass of NGC 6302's 
progenitor would lie between 4 M and 5 M Q to account 
for the transition in C/O ratio. This mass range is also 
confirmed by the observed gas phase abundances, and 
is in independent agreement with the comparison of the 
central star's temperature and the post-AGB tracks of 
Blocker (1995). Figure || shows the synthetic evolution 
of a 4.45 Mq model that matches the observational re- 
quirements. It can be seen that the star goes through 
a C transition over the last 2000 years. The final C/O 
ratio is as high as that found in circumstellar environ- 
ments which show the UIR bands (and not SiC, e.g. 
Barlow 1983, Roche 1989). 



* The lack of a continuous transition would preclude the ex- 
istence of type Ila PNe, with log(N/0)>- 0.6 (as defined by 
Faundez-Abans & Maciel 1987), which in Groenewegen &; de Jong 
(1993) are mainly obtained at high initial mctallicity, but are 
nonetheless observed to exist in the LMC (Leisy & Dennefeld 
1996). A limiting core mass of 0.8 Mq for HBB also constrains N 
enrichment to occur for stars close to the 2 nd dredge-up critical 
mass (Becker & Iben 1979), which does not reproduce the popu- 
lation of C-rich Type I PNe, except for progenitor masses greater 
than 6 M . 
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Figure 7: The upper plot shows the variation of the 
time spent in the C star phase at the end of the AGB 
as a function of initial mass. The lower two plots show 
the C/O and N/O ratios in the ejecta averaged over 
the last 10000 years of evolution, and the horizontal 
dashed lines correspond to gas phase measurements of 
the corresponding ratios (see text). 
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Figure 8: Schematic evolution of a 4.45 M Q star at the 
end of the AGB. The dotted line is the mass loss rate 
(left hand scale) and the solid line is the surface C/O 
ratio (right hand scale). The last thermal pulse changes 
the C/O chemical balance of the star's envelope, ap- 
proximately 2000 years before envelope ejection. 



Table 3: FWHM line- widths AV in kms -1 , together 
with the ionization potential (IP, eV) required to pro- 
duce the observed ion, and the wavelength A in fim. 



Br7 


[Mgvm] 


[Si vi] 


[Arvi] 


[Alv] 


A [/mi] 2.166 


3.028 


1.963 


4.528 


2.883 


IP [eV] 13.6 


225 


167 


75.0 


154 


NGC 6302 










AV | 31(5) 


18(2) 


22(5) 






AV_L 41(5) 


24(2) 




27(5) 




NGC 6537 










AV | 39(5) 




14(10) 




30(20) 



4 KINEMATIC AL STRUCTURE OF THE 
CORES OF NGC 6302 AND NGC 6537 

4.1 Velocity profiles from the echelle spectra 
and the observed ionization stage 
stratification 

The velocity profiles along the echelle slit are shown in 
Figure | for NGC 6302. A radically different velocity 
structure can be seen when comparing the coronal ions 
to Bi'7 for the _L slit orientation, along the bipolar axis. 
The ions show an outflow type of profile, whereas Br7 
displays the characteristic X-shaped profile obtained for 
optical recombination lines, and for [Nil] by Elliot & 
Meaburn (1977). This difference can be checked at offset 
-8", where the double gaussian fit to the Br7 line pro- 
file has clearly separated centroids, whereas [Mgvm] or 
[Arvi] are centrally peaked (the maxima at each offset 
are shown in crosses). The || slit orientation shows nar- 
row and centrally peaked velocity profiles, in particular 
for [Mgvm]. Br7 is less extended in the || direction, 
consistent with the structure expected from a bipolar 
flow in the _L direction. 

For NGC 6537, velocity profiles with the || slit orien- 
tation only are available. Figure [n^ shows that they have 
a similar structure to the || slit orientation in NGC 6302. 
The lines are also quite narrow. 

Table | lists the widths of the lines (FWHM, the in- 
strumental line widths were subtracted in quadrature) 
obtained by fitting the spectra extracted from the col- 
lapsed slits. The uncertainties stem from the instrumen- 
tal line- widths. The line- widths from the arc lamp spec- 
tra for the two [Mg vm] and the [Al v] echelle frames 
are consistent with AV=26±1 kms -1 over the central 
third of the array, although the [Ar vi] arc frame had 
a 20 kms -1 wide line in the centre. But it is clear that 
the lines are very narrow indeed, and might even not be 
resolved in the case of [Sivi] in NGC 6537. 

The velocity profiles give evidence for stratification 
in NGC 6302, as can be seen in Figure [H]. The con- 
tour at half maximum corresponding to [Mg vm] (dot- 
ted line), at IP= 224.95 eV, is less extended than those 
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Figure 9: Results of our echelle observations of NGC 6302. The two slit positions are referred to as || for the direction 
along the waist of the nebula, and _L along the bipolar axis. From left to right are shown Br7, [Mgvni] 3.028^m and 
[Arvi] 4.528^m. The y-axis is the offset along the slit in arcsec, and the x-axis shows the velocity in the object's rest 
frame (km s _1 ). The contour levels are at the maximum flux density times powers of 1/2. The crosses correspond 
to the maxima for each row. East is towards positive offsets in the frames with the _L slit orientation. 
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Figure 10: Results of our echelle observations of NGC 6537. From left to right arc shown Br7, [Si vi] 1.963/im and 
[Alv] 2.883/im. The slit was oriented along the waist of the nebula. The y-axis is the offset along the slit in arcsec, 
and the x-axis shows the velocity in the object's rest frame (km s _1 ). The contour levels are at powers of (1/2) times 
the maximum flux density for each frame. 
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Figure 11: Spatial and velocity extent of [Mgvm] (dot- 
ted line), [Sivi] and [Arvi]. The contours are taken at 
half maximum a) for the slit oriented along the dusty 
lane, b) for the slit oriented along the bipolar axis. The 
registration of the contours is uncertain to 0.1 — 0.2". 
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Figure 12: The ionization structure of a uniformly dense 
nebula as computed by CLOUDY, with NGC 6302's pa- 
rameters. 

of [Arvi], IP= 75.04 eV, and [Sivi], IP=166.77eV. The 
nebular emission in each frame did not necessarily peak 
on the same row, which affects the registration of the 
contours between different observations. The relative 
positioning of the contours is also affected by the limited 
spatial resolution of about 1.5" along the slit, and by 
the velocity resolution of ~10 km s _1 per pixel along the 
dispersion axis. It can be noted that the lower ioniza- 
tion species are more spatially extended than [Mg vm] , 
which matches the structure of a photoionizcd nebula, 
shown in Figure [j^. The velocity width of the lines also 
show stratification, implying that material is acceler- 
ated outwards from the nucleus. 



emerges from the optical studies is typically a bipolar 
cavity surrounded by an expanding H n region, which in 
turn is confined by some higher density medium (some- 
times detected in Hi or CO). The observational results 
are usually compared with the interacting stellar winds 
model (ISW, Kwok et al. 1978), with fairly good agree- 
ment. For instance, Corradi & Schwarz (1993) studied 
NGC 6537 and Hb 5 in [N n] and compared their results 
with the ISW model of Icke et al. (1989); NGC 6537 has 
also been studied in [Nil] by Cuesta et al. (1995) and 
its kinematical structure was consistent with refraction 
of a fast wind on the walls confining an evacuated cav- 
ity (model for Herbig-Haro objects from Canto 1980); 
Bryce et al. (1996) established the velocity structure of 
NGC 650-1 using [Nil] and [Oin], and the results were 
in reasonable agreement with the GISW models pre- 
sented in Frank and Mellcma (1994). A common feature 
in the models quoted previously is that a hot evacuated 
cavity (N c ~ 10 cm -3 ) occupies a good portion of the 
PN's volume, whose walls are expanding with velocities 
of order 50 kms -1 . 

In the case of NGC 6302, the above picture would 
seem to be confirmed by the CGS4+echelle observations 
of Bra. Fig. [9] shows that the Bra profiles double-peak 
further out than 8" from the nucleus, which is an in- 
dication for a possible decrease in proton density along 
the bipolar axis. This could also be extrapolated from 
the optical echelle measurements of recombination lines 
(Elliott & Meaburn 1977), and from the H/3 image of 
Ashley (1990), under the assumption of axial symmetry. 

But the use of coronal ions as tracers of the ve- 
locity field produces quite a different picture. The ob- 
served ionization stage stratification is consistent with 
a CLOUDY model for a nebula at constant density (Fig- 
ure |l^). We also have no indications for a torus type 
of structure rather than a disk around the waist of the 
nebula, at least at about 1" resolution, as can be seen 
from the frames along the dusty lane. The same situa- 
tion seems to apply to NGC 6537 (Figure |ic|). Addition- 
ally, the single peaked profiles of the coronal lines out to 
8" offsets along the bipolar axis (in NGC 6302, Fig. |) 
show that the corresponding physical conditions cannot 
be that of a rarefied plasma. The detection of coronal 
line emission requires N p > 10 3 cm -3 out to 10" offsets, 
and even higher if the decrease in ionic abundance is 
taken into accountfj]. 

Arguments against the existence of a sizeable evac- 
uated cavity may also be extrapolated from the VLA 
6 cm map of NGC 6302's central 20" by Gomez et al. 



4.2 The absence of a wind blown cavity 

Most kinematical studies of bipolar nebulae have been 
based on optical echelle spectroscopy of [Nil] 6584 A, 
which does not probe the high excitation regions and 
suffers from extinction. The general structure that 



T Using detailed balance it is possible to determine the range in 
temperature and density where the predicted ratio of the flux at 
nebular conditions (N e = 18000 cm -3 and T c = 20000 K, Section 
0) to the flux at (N e ,T) would put it below the noise level. A drop 
in [Mgvm] flux by a factor 1000 corresponds to a drop in density 
from nebular conditions (Section ^) to N p > 8 X 10 2 cm~ 3 . The 
limiting case N p = 8 X 10 2 cm -3 corresponds to T e = 10 4 K. 
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(1993). The 6cm continuum traces the emission mea- 
sure across the central regions of the nebula (provided 
T c is fairly uniform). Under the assumption of axial 
symmetry, the possibility of an extreme decrease in pro- 
ton density from nebular conditions to ~10 cm~ 3 along 
the bipolar axis can be discarded (note that for a given 
6 cm flux density an increase in T requires an increased 
EM) . The CGS4+echelle results reported here confirm 
the filled-in picture evident from the VLA map, this 
time using tracers of highly excited gas. Also early ob- 
servations of PNe by Wilson (1950) showed that the cen- 
tral 'cavities' in the long slit spectra of low-ionization 
species (such as [O n] and [N n] ) are filled-in, spatially 
and spectrally (at ~10 kms" 1 resolution), by highcr- 
ionization species such as [Nev]. Wilson (1950) thus re- 
ported a stratification in ionization potential, including 
the intermediate ionization species [O in] and [Ne in] . 

NGC 6302 has previously been studied in [Nev] by 
Meaburn & Walsh (1980). They concentrated on the de- 
tection of a fast wind in the broad wings of [Ne v] , which 
are absent in their [N n] profiles. This discovery has been 
interpreted by Barral et al. (1982) in terms of the Canto 
(1980) model for H-H objects. In this model, flows are 
generated as the stellar wind refracts on a dense disk 
surrounding the star and is accelerated along the walls 
of the bipolar nebula. The shape of the walls is deter- 
mined by equating the ram pressure of the wind to that 
of the surrounding medium. The conclusion of Barral 
et al. (1982) within this model is that the dominant ex- 
citation mechanism is still radiative ionization. In this 
case, an altogether different model for photoionization 
should be considered, in which the nebula has in fact 
a shell structure, the hot bubble consisting only of the 
stellar wind itself. But the spectra of [Nev] published 
by Meaburn & Walsh do not show a double peaked 
structure, at a resolution of about 10 kms -1 . Although 
Meaburn & Walsh (1980) interpreted the geometry of 
NGC 6302 as an evacuated cavity, with [Ne v] emission 
originating from the walls, their results do not show 
the corresponding velocity profile. We believe a filled- 
in structure is not inconsistent with observing different 
velocity components, provided the nebula has a clumpy 
structure. 

This Section points towards the absence of a wind- 
blown cavity filled with a hot rarefied plasma, in the 
case of the highest excitation PN known. We suggest 
that the dumpiness of the nebula precludes the for- 
mation of a wind-blown cavity. This point should be 
further investigated to see under which conditions, if 
any, bipolar PNe would show signs of a hot bubble. 
Similar conclusions have been reached by O'Dell (1998) 
and Henry et al. (1999) in the case of the Helix Nebula 
(NGC 7293), namely that the ionized phase of the neb- 
ula is a disk rather than a ring, filled by highly ionized 
gas. The current ISW dynamical models that account 
for the shaping of PNe have also been used to explain 
a variety of astrophysical objects, and even the collima- 



tion of jets in AGN (Icke et al. 1992) and young stellar 
objects (Mellema & Frank 1997). But the ISW model 
was built from and checked against optical studies of 
relatively low ionization species in PNe, and we have 
shown here that the observational information should 
be extended by studies of the cores of other nebulae in 
higher-excitation species. 

4.3 What are the physical conditions in the 
fast wind? 

The signal to noise ratio of our echelle observations of 
NGC 6302 in [Ar vi] and [Mg vm] is higher than that of 
[Nev] 3426 A presented by Meaburn & Walsh. Also, 
Nev has an ionization energy intermediate between 
that of Arvi and Mgvni, while the IR coronal lines 
and the optical [Ne v] line have similar critical densi- 
ties (at 10 4 K, the IRON project collision strengths give 
A crit ([Arvi])~ 10 6 cm- 3 , A crit ([Ne v])~ 5xl0 6 c m - 3 , 
-Nci-it ([Mgvin])~ 10 7 cm -3 ). So why are no extended 
wings observed in the IR coronal lines? 

Without mass loading, the fast stellar wind should 
vary in density with angular distance in arcsec, a, from 
the central star approximately as 10/a 2 cm -3 /" 2 (for 
dM/dt = IQ- S Moyr" 1 and V lw = 1000 kms" 1 ). This 
would not be detected in IR emission because of insuf- 
ficient collisional excitation, leading to a very low ex- 
cited ion column density. But the physical conditions 
where the broad wings originate could be intermediate 
between the fast stellar wind and nebular conditions. 
Let us consider the implications of assuming collisional 
excitation for the fast wind observed by Meaburn & 
Walsh (1980). For the 3p 2 - 1 D 2 [Nev] 3426A transition, 
the line intensity can be written as 

h-j = -^Jedsb A i:j hv n up N a , (2) 

where b is the branching ratio, e the filling factor, 
n up is the fraction of ions in the upper level, and N a 
is the number density of ions. We have given a number 
of justifications for the fast wind coexisting with the 
denser clumps, which occupy a significant fraction of 
the volume (e ~ 0.5), so the fast wind's filling factor 
should be similar. The ratio of the intensity from the 
narrow component of the line to that from the fast wind 
will be given approximately by 

-^narrow ^ (^up-^p)narrow 
-Swings (^up-Mp)wings 

as long as the ionic abundance does not change signifi- 
cantly between the two components. Meaburn & Walsh 
(1980) give 

,-wings /o-, riiarrow ( a\ 

Vev] ~ ( 7 j [Nev]" ^ 

If the ionic abundances were substantially different in 
the fast wind, such an important contribution to the 
total line flux would put into question photoionization 
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Figure 13: Estimation of the ratio of the broad com- 
ponent to the narrow component as a function of 
log(T/10 4 K) (ordinate) and electron number density. 
The transitions correspond to [Ne v] 3426A and [Ar Vl] 
4.528^tm, and [Mgvm] 3.0284^im. The thick line limits 
the locus of {log(T/10 4 K), N e } that can account for the 
broad wings observed in [Ne v] 3426A . 

models. A two-phase nebular gas should be considered. 
Our IR observations do not show the fast wind, so the 
photoionization models presented here should not be 
affected since we use no optical lines. A detailed pho- 
toionization model including a range of optical and in- 
frared lines could constrain possible variations in ionic 
abundances among the two phases. Nonetheless, under 
the assumption that ionic abundances do not change be- 
tween clump and fast wind, the range of physical condi- 
tions in the fast wind that gives a [Ne v] line ratio of 1 /3 
is shown in Figure [r], together with the corresponding 
ratio for the case of [Ar vi] and [Mg vm] . The noise lev- 
els reported in Figure |^ yield the following upper limits 
on the broad wings fluxes in [Mgvm] and [Arvi], as- 
suming the broad wings flux is uniformly distributed 
over 800 kins" 1 , 



jwings ^ 

[Mg vm] 



fl HO rnarrow 

U ' U ^ 1 [Mgvm] 



7-wings ^ r narrow 

Vrvi] ~ U '^ 7 [Arvi] 



(5) 
(6) 



Using Figure it can be concluded that broad wings 
should also be apparent in [ArVl], and especially in 
[Mg vm] , for the acceptable values of densities and tem- 
peratures. It seems that collisional excitation cannot 
account for the broad wings observed in the optical, 
at least in these simple terms. A more thorough study 
including a broader range of lines is required. 



5 CONCLUSIONS 

The high excitation PNe NGC 6302 and NGC 6537 
have been studied with medium and high resolution in- 
frared spectroscopy. We provide a new IR estimate for 
the heliocentric radial velocity of NGC 6302, finding a 



4.1 ± 2.7 kms -1 difference from optical estimates, pos- 
sibly due to extinction within the ionized region. The 
coronal line wavelengths we derive are in agreement 
with those of Reconditi & Oliva (1993), but not with 
those of Feuchtgruber et al. (1997). 



Our medium resolution spectra show 
[Si ix] 3.934/im in emission, which is used with 
other IR coronal ions to confirm the nebular excitation 
mechanism to be photoionization, corresponding to a 
blackbody central star temperature of about 250 000 K 
for NGC 6302. The only anomalous elemental abun- 
dance derived from our NGC 6302 spectra is that of 
aluminium, found to be depleted relative to solar by a 
factor of more than one hundred, presumably due to 
being locked in grains. The strong Al depletion implies 
that some dust is well mixed with the highly ionized 
gas where the radiation field is hard enough to produce 
Alvi, at IP=154 eV. A similar Al depletion is inferred 
for NGC 6537. For the case of single star evolution, 
the revised central star temperature for NGC 6302 is 
consistent with post-AGB model tracks (Blocker 1995) 
for ZAMS masses between 4 and 5 M . NGC 6302's 
N/O and C/O gas phase abundance ratios are also 
consistent with a progenitor main sequence mass < 5 
M Q , from an adaptation of the synthetic AGB models 
of Gronewegen & deJong (1993). 



Our IR echelle observations reveal that the infrared 
coronal lines are very narrow and strongly peaked at 
the centre, with no evidence of broad wings from fast 
winds. The discrepancy with the ~1000 kms" 1 broad 
wings previously observed for [Ne v] 3426 A cannot be 
explained by collisional excitation in a photoionized 
gas, at least in a simplified treatment. Further obser- 
vations are required to confirm the [Ne v] broad wings, 
and to build detailed models of the interaction be- 
tween the fast stellar wind and the nebula. The nar- 
rowness of the line profiles also makes NGC 6302 and 
NGC 6537 ideal objects for the study of fine-structure 
in coronal ions. The high resolution observations show 
that the cores are stratified in ionization stages, con- 
sistent with a photoionized nebula of uniform density. 
At 8" offset from the centre of NGC 6302, Br 7 shows 
double-peaked profiles, in agreement with previous op- 
tical studies based on recombination lines and [Nil]. 
In contrast, the line profiles from the coronal ions are 
single-peaked at the same offset. Thus our observations 
show no signs of a wind-blown bipolar cavity. The neb- 
ula, although clumpy enough to accommodate coexist- 
ing winds at very different velocities, seems to have a 
fairly uniform density as close as 1" to the nucleus. We 
suggest that the dumpiness of the nebula precludes the 
formation of a wind-blown cavity. 
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APPENDIX A: FLUXES AND ACCURATE 
WAVELENGTHS FOR FINE STRUCTURE 
LINES, AND THE HELIOCENTRIC 
VELOCITIES OF NGC 6302 AND NGC 6537 

A number of fine structure lines can be found in the 
spectrum of NGC 6302 and NGC 6537, whose accurate 
wavelengths represent valuable information and can be 
compared with previous measurements. The low reso- 
lution spectra were calibrated using nebular Hi lines 
when present near the lines of interest, otherwise using 
a quadratic fit to Ar or Kr arc lamp spectra (deviations 
from a linear law were found to be negligible). The cal- 
ibration of the echelle observations is described below. 
We found the heliocentric velocity of NGC 6302 to be 
-34.8±1 kms -1 , using Hi 5-10 at 3.0384/im in the two 
[Mg vni] frames. This corresponds to -27.3±1 kms -1 in 
the local standard of rest (kinematic LSR), which agrees 
closely with the H76a measurement of -29±2 kms -1 by 
Rodriguez et al. (1985, their Fig. 2 shows that the ob- 
served peak is red shifted with respect to their Gaussian 
fit by about 5 kms -1 ). But the heliocentric velocity 
of NGC 6302 obtained from optical studies (compila- 
tion by Durand et al. 1998, based on Acker 1975) has 
been reported to be -38.9±2.5 kms -1 . On the other 
hand we measure a heliocentric velocity for NGC 6537 
of-17.8±3 kms -1 , using Hi 5-11 at 2.872213/mi on the 
[A1V] frame, which is in agreement with the value listed 
in Durand et al. (1998) of -16.9±1.9 kms" 1 . It should 
be stressed that the long wavelength side of Hi 5-11 
was very noisy, the heliocentric velocity we report for 
NGC 6537 assumes that the width of the line is the 
same as for Br7, but any error would result in a higher 
(i.e. red shifted) velocity. 

By adopting the estimated systemic velocities from 
the echelle spectra, we can achieve consistent wave- 
length determinations for the extinction insensitive in- 
frared coronal lines. Considering the quoted uncertain- 
ties in the case of NGC 6302, there could be a slight 
discrepancy with the optical measurements, which is 
most likely due to internal extinction within the ion- 
ized regions of the nebula. Huggins & Healy (1989) 
report a CO(2-l) centroid for NGC 6302 of V, sr = 
40 kms -1 , blue shifted in excess of 10 kms -1 from 
the measurements reported here. The CO(2-l) profile 
is double peaked, and the approaching side has about 
twice the peak brightness temperature as the reced- 
ing side (we conducted SEST observations of 13 CO(2-l) 
which show the same profile, and that shall be presented 
in a forthcoming article), hinting at possible radiative 
transfer effects. However, we do not base any of the 
following on the mm-line centroids, because the line of 
sight towards NGC 6302 crosses cold galactic molecular 
emission. 



Table Al 



gives a summary of our measurements, 
where wavelengths are given in air and in the nebu- 
lar rest frames (the refractive index of standard air as 



a function of wavelength was taken from Edlen 1953). 
The uncertainty in the calculated wavelengths is dom- 
inated by the accuracy of the adopted heliocentric ve- 
locities, in the best cases. Below we provide some notes 
and compare our measurements with those found in the 
literature, where available. 

[Mgvili]: We have two echelle measurements cor- 
responding to the two slit orientations. A comparison 
of the standard star spectra with a high-resolution at- 
mospheric transmission model (based on the HITRAN 
database, Rothman et al. 1992) gave an initial quadratic 
fit to the dispersion law, which we used to identify Ar- 
gon lines from the arc lamp frames. The arc lamp cal- 
ibration agrees very well with the atmosphere calibra- 
tion, although with shifts of order 10 kms -1 towards the 
edges of the array, which probably stem from the stan- 
dard star not being precisely centred within the 1" slit. 
We used the arc lamp to calibrate the [Mgvin] frame. 
Many lines were identified as leaks from neighbouring 
orders, i.e. the Argon arc lamp spectra in order 36 also 
showed lines from orders 37 and 38, implying that the 
CVF blocked wavelengths outside of about 0.1/im from 
the CVF wavelength of 1.515 /im. The calibrated disper- 
sion axis is strongly non-linear, with deviations from a 
linear law of up to 20 kms -1 , and offset from the es- 
timated scale by -300 kms -1 to -50 kms -1 from side 
to side of the array. We estimate the uncertainty in the 
[Mgvin] wavelength to be 2 kms -1 , dominated by the 
relative velocity of NGC 6302. Hi 5-10 is on the same 
frame, and allows an accurate estimate of the doppler 
shift from NGC 6302. The low resolution spectrum es- 
timate is accurate to better than one sample element 
(i.e. R=1500x3), but the secondary line is blended. 

[ArVi]: The echelle wavelength calibration was 
based on a xenon lamp, which showed one arc line only, 
at 2.2618283/im (Outred 1978) in second order, very 
close to the centre of the array. Because of the mod- 
ulation due to the CVF fringes it was impossible to 
use the standard star spectrum to calibrate the wave- 
length dispersion. We used the quadratic dispersion law 
for [Mgvin] 3.03 /im, shifted to match the xenon arc 
line. The uncertainty should be about 10 kms -1 . The 
low resolution spectrum did not provide a satisfactory 
measurement because the line was strongly affected by 
random bad pixels, with an accuracy of about one res- 
olution element. 

[Al v]: Here the echelle spectrum was of NGC 6537, 
rather than NGC 6302. We calibrated the dispersion 
axis by the same procedure as for [Mgvin], and con- 
firm that the CVF does not block wavelengths within 
0.1 /mi from the CVF wavelength. Hi 5-11 2.872213/mi 
was also present in the same frame, although it is noisy 
this line allows an estimate of the systemic velocity of 
NGC 6537. No measurements are provided for the low 
resolution spectrum because the line is blended with H I 
5-11, which represents about 75 per cent of the flux. 

[Sivi]: we used the airglow OH line at 1.95878 /im 
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Table Al: Rest wavelengths (in air) for emission lines from NGC 6302, or from NGC 6537 when specified. The 
numbers in brackets are the uncertainties in the last decimal places. Fluxes are given for the whole slit, without 
correction for extinction or blends and with a typical uncertainty of 20 per cent. The slit orientation for each line 
can be inferred from the tabulated fluxes. 

Ion i < j A [/im] A [/xm] A [/xm] 6 Flux [W m" 2 ] 

R=1500 R=20000 (reference) || 3" slit || 1" slit _L1" slit 



Si vi 


2 P 3 /2- 


■ 2 Pl/2 




1.96287(10) 


1.96287(10) 


(3) 






1.6x10" 


-16 




Si vi c 


2 P 3 /2- 


- 2 Pl/2 




1.96311(20) 


1.96287(10) 


(3) 






1.0x10" 


-17 




Alv c 


2 P 3 /2- 


- 2 Pl/2 




2.882975(19) 


2.90440(37) 


(1) 






9.6x10" 


-18 




Mg vin 


-tl/2" 


2p 

' M/2 


3.02788(67) 


3.027661(20) 


3.02713(25) 


(3) 


2.9x10" 


-15 


1.7x10" 


-15 




(second) a 




3.027033(20) 










1.9x10" 


-16 




Mg VIII 


2-p 
-tl/2" 


2p 

' ^3/2 




3.027648(20) 


3.02713(25) 


(3) 










1.4xl0- 15 


(second) a 




3.027003(20) 














7.3xl0"" 17 


Kvn 


2 Pl/2- 


- 2 P 3 /2 


3.1899(21) 




3.18966(15) 


(1) 


4.2x10" 


-16 








Alvi 


3 P 2 - 


■ 3 Pl 


3.65952(81) 




3.661(14) 


(2) 


2.1x10" 


-16 








Si ix 


3 Po- 


- 3 Pl 


3.9346(26) 




3.935 


(4) 


3.5x10" 


-17 








Arvi 


2 Pl/2- 


- 2 P 3 /2 


4.5265(30) 


4.52799(15) 


4.52829(31) 


(1) 


8.7x10" 


-14 






2.5xl0" 14 


Na vii 


2 Pl/2- 


- 2 P 3 /2 


4.6908(31) 




4.68344(33) 


(1) 


7.4x10" 


-15 









a Unidentified secondary line in the [Mgvm] 3.028/im frame. 

b Previous measurements from 1), Feuchtgruber et al. (1997, corrected to values in air); 2, Greenhouse et al. (1993); 
3), Reconditi & Oliva (1993); 4), Oliva et al. (1994). 
c Measured from NGC 6537. 



(1.95932 fim in vacuum from the UKIRT worldwide 
webpages). The line is very faint and was discernible 
across the array only after co-adding all the object 
frames, and there is an uncertainty of about 10 kms 
in the position of the line. The standard star spectra 
seem to be affected by local absorption related to the 
instrument, which is also apparent in the nebular con- 
tinuum emission from both PNe. The solar spectrum is 
at odds with the standard's, so unfortunately the at- 
mosphere cannot be used for calibration. We adopted 
the dispersion scale from the [Mgvm] frame shifted 
to match the OH line. An uncertainty of no less than 
10 kms" 1 is expected in the dispersion scale, and the 
[Sivi] wavelength is thus accurate to about 15 kms" 1 . 

[Alvi]: Only a low resolution wavelength determi- 
nation is available, but it should be reliable because 
that spectrum showed many other lines, which all match 
the reference wavelengths. However the best accuracy 
achievable is probably not better than one sample el- 
ement, or 67 kms" 1 (R = 1500 oversampled 3 times), 
because the line is noisy and lies near a zone of strong 
atmospheric absorption. 

[Siix]: this line was calibrated using nearby Hi 
Hul5, but it is faint and the accuracy is only one resolu- 
tion element. The resulting wavelength is in good agree- 
ment with the measurement from the coronal spectrum 
of the Circinus galaxy by Oliva et al. (1994). This comes 
as a confirmation for their adopted radial velocity since 
in Circinus the narrow line region is blue shifted relative 



to the systemic velocity, and [Siix] is about 300 kms" 1 
broad. 

Finally, [Kvn] was calibrated with respect to an 
Argon arc lamp and [Na vii] with respect to the nearby 
Hen 6-7 line, although both were quite noisy and so 
the wavelengths are estimated to be accurate to only 
one resolution element. 

The low resolution and the echelle measurements 
are consistent, and agree in general with other published 
results, except for the ISO measurements of Feuchtgru- 
ber et al. (1997). The agreement with Reconditi and 
Oliva (1993) is improved when using the optical helio- 
centric velocity of NGC 6302 discussed above. 

The fluxes for emission lines in the CGS3 spectra 
are listed in Table |A^, together with the ISO fluxes from 
Pottasch & Beintema (1997) and Pottasch et al. (1996). 
In order to compare the CGS3 and ISO fluxes, correc- 
tion factors need to be applied to account for the differ- 
ent beam sizes, which can be estimated using the spa- 
tial information from the low resolution CGS4 spectra. 
In the case of NGC 6302, at a 3" pixel size the average 
spatial full width at half maximum (FWHM) of HI lines 
along the slit was 2.8 pixels, and all coronal ions had ap- 
proximately the same width of 2 pixels. For NGC 6537, 
the HI FWHM was 1.7 pixels, and the coronal ions cov- 
ered approximately 1.3 pixels along the slit. These re- 
sults suggest that, for the purpose of scaling fluxes to 
the central 3", all ions can be approximated to have the 
same spatial extent. At the higher spatial resolution of 
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Table A2: Fluxes of emission lines in the 8-13/xm spec- 
tra of NGC 6302 and NGC 6537 (in 10~ 14 Wm~ 2 ), and 
comparison with ISO measurements. 







IP 


NGC 6302 


NGC 6537 






[eV] 


CGS3 


ISO 


CGS3 


[Ne vi 


7.64 /im 


126.1 


81.9 


63.0 


13.5 


[Ar v] 


7.90 fim 


59.8 


1.18 




0.86 


[Navi 


8.62 jurn 


138.4 


0.56 


1.1 


0.23 


[Ar in 


8.99 ^m 


27.6 


1.76 




0.80 


[Siv] 


10.5 /urn 


34.8 


4.77 


8.5 a 


6.77 


[Neil] 


12.8 /xm 


21.6 


3.23 


14.7 a 


1.14 


[Ar v] 


13.1 fj,m 


59.8 


0.56 


2.9 


0.61 


from 


Pottasch 


& Beintema (1997), the other ISO 



fluxes are from Pottasch et al. (1996). 

the echelle observations (see Figure |ll|) the FWHM for 
[Mgvm], [Sivi] and [Arvi] are about 3", 4" and 5", 
and ions with ionization potentials ^100 eV would have 
a spatial extent intermediate between Sivi and Arvi. 
The fluxes from large beam observations should thus 
be about 2-3 times higher than from the central 3". 
This seems reasonable for [Mgvm] 3.028 /im, listed in 



Table Al, for which Pottasch et al. (1996) give a flux 
of 3-6xl0~ 15 Wm- 2 . For NGC 6302, fluxes through 
the circular 4" aperture of CGS3 should thus be 2- 
3 times less than the fluxes for the whole nebula, as- 
suming the emission lines have uniform surface bright- 



ness over the central square 5-6". From Table A2, a 
correction factor of ~ 2.5 matches [Navi] 8.62 /im and 
[Siv] 10.52 fim. We attribute the larger correction fac- 
tor required for [Ne n] 12.8 /im to a larger spatial extent, 
but the [Ne vi] 7.64 /Ltm fluxes are discrepant. [Ne vi] lies 
on the edge of the atmospheric window, and is likely to 
be affected by unresolved atmospheric transmission fea- 
tures. 



This paper has been produced using the Royal Astro- 
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